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Abstract 

The theoretical interest in the longitudinal diffractive structure function F[' is 
briefly motivated and possible measurement methods are surveyed. A simula- 
1/^ ■ tion based on realistic scenarios with a reduced proton beam energy at HERA- 

O I II using the HI apparatus shows that measurements are possible with up to 4a 

^ I significance, limited by systematic errors. 

> 

O ■ 1 Introduction 

I In order to understand inclusive diffraction fully, it is necessary to separate out the contributions from 

^ ' transversely and longitudinally polarised exchange photons. Here, the formalism of [1] is adopted, where 

by analogy with inclusive scattering and neglecting weak interactions, a reduced cross section aj? is 
defined,^ related to the experimentally measured cross section by 

> ' 

and y+ = 1 + (1 — y)^. The structure function F[', is closely related to the longitudinal photon con- 



o 



tribution, whereas the more familiar F2 contains information on the sum of transverse and longitudinal 
photon contributions. 



^ ■ 9 n 

D , It is generally understood [2] that at high /? and low-to-moderate Q , receives a significant, per- 

Q_i' haps dominant, higher twist contribution due to longitudinally polarised photons. Definite predictions [3] 

1^ ■ exist for this contribution, obtained by assuming 2-gluon exchange, with a similar phenomenology to that 

' successfully applied to vector meson cross sections at HERA. The dominant role played by gluons in the 

diffractive parton densities [1] implies that the leading twist F£ must also be relatively large. Assuming 
^ ■ the validity of QCD hard scattering coUinear factorisation [4], this gluon dominance results in a leading 

^ ', twist F[' which is approximately proportional to the diffractive gluon density. A measurement of F[' 

to even modest precision would provide a very powerful independent tool to verify our understanding 
of the underlying dynamics and to test the gluon density extracted indirectly in QCD fits from the scal- 
ing violations of ^2^. This is particularly important at the lowest x values, where direct information on 
the gluon density cannot be obtained from jet or D* data due to kinematic limitations and where novel 
effects such as parton saturation or non-DGLAP dynamics are most likely to become important. 

Several different methods have been proposed to extract information on F^. It is possible in 
principle to follow the procedure adopted by HI in the inclusive case [5,6], exploiting the decrease in 
at large y relative to expectations for alone (see equation 1). This method may yield significant 
results if sufficient precision and y range can be achieved [7], though assumptions are required on the xp 
dependence of ^2^, which is currently not well constrained by theory. An alternative method, exploiting 
the azimuthal decorrelation between the proton and electron scattering planes caused by interference 
between the transverse and longitudinal photon contributions [8], has already been used with the scattered 
proton measured in the ZEUS LPS [9]. However, due to the relatively poor statistical precision achievable 
with Roman pots at HERA-I, the current results are consistent with zero. If the potential of the HI VFPS 
is fully realised, this method may yet yield significant results in the HERA-II data [10]. However, if 
the necessary data are taken, the most promising possibility is to extract F[' by comparing data at the 

'it is assumed here that all results are integrated over t. The superscript (3) usually included for F,^'^^^ and other quantities 
is dropped for convenience. 



same Q^, (3 and xjp, but from different centre of mass energies y/s and hence from different y values. 
The longitudinal structure function can then be extracted directly and model-independently from the 
measured data using equation 1. In this contribution, one possible scenario is investigated, based on 
modified beam energies and luminosities which are currently under discussion as a possible part of the 
HERA-II programme. 

2 Simulated Measurement 

Given the need to obtain a large integrated luminosity at the highest possible beam energy for the re- 
mainder of the HERA programme and the fixed end-point in mid 2007, it is hkely that only a relatively 
small amount of data can be taken with reduced beam energies. A possible scenario is investigated here 
in which 10 pb^^ are taken at just one reduced proton beam energy of Ep = 400 GeV, the electron 
beam energy being unchanged at 27.5 GeV. Since the maximum achievable instantaneous luminosity at 
HERA scales like the proton beam energy squared [1 1], this data sample could be obtained in around 2-3 
months at the current level of HERA performance. It is assumed that a larger data volume of 100 pb^^ is 
available at Ep = 920 GeV, which allows for downscaling of high rate low inclusive triggers.^ The 
results presented here can be used to infer those from other scenarios given that the statistical uncertainty 
scales like a^"^^^ /\fC^ + u^^"^^ /\fC^, where a^^^ and are the reduced cross section and the 
luminosity at a proton beam energy of Ep, respectively. 

The longitudinal structure function can be extracted from the data at the two beam energies using 



where and Y_^'^ denote y and i+ at a beam energy Ep. It is clear from equation 2 that the best 
sensitivity to requires the maximum difference between the reduced cross sections at the two beam 
energies, which (equation 1) implies the maximum possible y at Ep = 400 GcV. By measuring scattered 
electrons with energies E'^ as low as 3 GeV [5], the HI collaboration has obtained data at y = 0.9. This is 
possible with the use of the SPACAL calorimeter in combination with a measurement of the electron track 
in either the backward silicon tracker (EST) or the central jet chamber (CJC). For HERA-II running, the 
corresponding available range of scattered electron polar angle is 155° < 9'^ < 173°, which is used in the 
current study.^ Three intervals in y are considered, corresponding at Ep = 400 GeV to 0.5 < y^Qo < 0.7, 
0.7 < j/400 < 0.8 and 0.8 < y4oo < 0.9. It is ensured that identical ranges in j3, xp and axe studied 
at Ep = 920 GeV by choosing the bin edges such that 7/920 = 2/400 " 400/920. Since the highest 
possible precision is required in this measurement, the restriction xjp < 0.02 is imposed, which leads 
to negligible acceptance losses with a typical cut on the forwardmost extent of the diffractive system 
Vmax < 3.3. The kinematic restrictions on E'^, 6'^ and xp lead to almost no change in the mean Q^, 
or /? ~ Q^/ (Q^ -I- M^) as either y or Ep are varied. In contrast, xp = Q^/{sy (3) varies approximately 
as 1/y. As is shown in figure 1, at the average (3 = 0.23, there is at least partial acceptance for all y bins 
in the range 7 < < 30 GeV^, which is chosen for this study, leading to an average value of close 
to 12 GeV^. 

The simulation is performed using the RAPGAP [13] Monte Carlo generator to extract the number 
of events per unit luminosity in each bin at each centre of mass energy. The values of F2' and F^, and 



hence and are obtained using an updated version of the preliminary HI 2002 NLO QCD 

fit[l]. 



^Alternative scenarios in which a smaller data volume at large Ep is taken in a short, dedicated run, could potentially lead 
to better controlled systematics at the expense of increased statistical errors. 

^One interesting alternative running scenario [12] is to obtain data at Ep = 920 GeV with the vertex shifted by 20 cm in 
the outgoing proton direction, which would allow measurements up to = 175°, giving a low acceptance range which 
closely matches that for the Ep = 400 GeV data at the normal vertex position. 
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Fig. 1: Illustration of tiie kinematic plane in and xp at proton energies of 920 GeV and 400 GeV, with fixed j3 = x/xp = 
0.23. The solid lines illustrate the experimental limits of 155° < 6'^ < 173°. The horizontal dashed lines illustrate the 
range used for the simulation. The diagonal dashed lines illustrate the binning in y, corresponding at Ep = 400 GeV toy = 0.9 
(leftmost Une), y = 0.8, y = 0.7 and y = 0.5 (rightmost line). 

The expected precision on is obtained by error propagation through equation 2. The system- 
atic uncertainties are estimated on the basis of previous experience with the HI detector [1,5]. At the 
large y values involved, the kinematic variables are most accurately reconstructed using the electron en- 
ergy and angle alone. The systematic uncertainties on the measurements of these quantities are assumed 
to be correlated between the two beam energies. With the use of the BST and CJC, the possible bias in 
the measurement of 9'^ is at the level of 0.2 mrad. The energy scale of the SPACAL calorimeter is known 
with a precision varying linearly from 2% at E'^ = 3 GeV to 0.2% at E'^ = 27.5 GeV. Other uncer- 
tainties which are correlated between the two beam energies arise from the photoproduction background 
subtraction (important at large y and assumed to be known with a precision of 25%) and the energy scale 
for the hadronic final state used in the reconstruction of and hence xjp (taken to be known to 4%, as 
currently). Sources of uncertainty which are assumed to be uncorrelated between the low and high Ep 
measurements are the luminosity measurement (taken to be ±1%), the trigger and electron track efficien- 
cies (±1% combined) and the acceptance corrections, obtained using RAPGAP (±2%). The combined 
uncorrelated error is thus 2.4%. Finally, a normalisation uncertainty of ±6% due to corrections for pro- 
ton dissociation contributions is taken to act simultaneously in the two measurements. Other sources 
of uncertainty currently considered in HI measurements of diffraction are negligible in the kinematic 
region studied here. 
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Table 1: Summary of the simulation at = 12 GeV and /3 = 0.23. The first three columns contain the y ranges used at Ep = 
400 GeV and Ep = 920 GeV and the xp values. The next two columns contain the values of the diffractive structure functions. 
These are followed by the uncorrelated ((5unc) and proton dissociation (Jnorm) uncertainties and the correlated systematics due to 
the electron energy (5E',) and angle (S9',) measurements, the hadronic energy scale (SAI^ ) and the photoproduction background 
(S'yp), all in percent. The last three columns summarise the systematic, statistical and total uncertainties. 
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Full details of the simulated uncertainties on the 
measurements are given in table 1. An illustration of the 
corresponding expected measurement, based on the Fj^ 
from the HI 2002 fit is shown in figure 2. The most pre- 
cise measurement is obtained at the highest y, where FJ^ 
would be determined to be unambiguously different from 
its maximum value of F2' and to be non-zero at the 4a 
level. Two further measurements are obtained at lower y 
values. The dominant errors arise from statistical uncertain- 
ties and from uncertainties which are uncorrelated between 
the two beam energies. Minimising the latter is a major ex- 
perimental challenge to be addressed in the coming years. 



0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.2 0.21 0.22 



Fig. 2: Illustration of the simulated result for Fj^ , 
showing the three data points with statistical (inner 
bars) and total (outer bars) errors. 



Only one possible scenario has been investigated 
here, leading to a highly encouraging result at relatively low 
/3, which would provide a very good test of the leading twist 
and thus of the gluon density extracted in QCD fits to 
^2^. It may also be possible to obtain results at high (3, 
giving information on the higher twist contributions in that 
region, for example by restricting the analysis to lower xp. 
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